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Abstract 

The generating function of Stieltjes-Carlitz polynomials is a solution of Heun's differ- 
ential equation and using this relation Carlitz was the first to get exact closed forms 
for some Heun functions. Similarly the associated Stieltjes-Carlitz polynomials lead 
to a new differential equation which we call associated Heun. Thanks to the link 
with orthogonal polynomials we are able to deduce two integral relations connecting 
associated Heun functions with different parameters and to exhibit the set of asso- 
ciated Heun functions which generalize Carlitz's. Part of these results were used by 
the author to derive the Stieltjes transform of the measure of orthogonality for the 
associated Stieltjes-Carlitz polynomials using asymptotic analysis; here we present 
a new derivation of this result. 



1 Introduction 



& ' The Stieltjes-Carlitz polynomials are othogonal polynomials with a three term recurrence 



relation 



with the rates 



(A n + Hn- x)F n (x) = fi n+1 F n+1 + A n _iF n _i, n > 0, 

F.i(x) = 0, F {x) = l, { " ! 



X n = k 2 (2n + 1) 2 , fi n = (2n) 2 

X n = (2n + 1) 2 , fi n = k 2 (2n) 2 

and < k 2 < 1. Their orthogonality measure, given for instance in || p. 194], was 
first derived by Stieltjes JL3| using continued fraction techniques. Later on Carlitz [0 



obtained a generating function for the F n (x) from which he was able to derive anew the 
orthogonality measure. The most striking fact, which does not seem to have been realized 
by Carlitz himself, is that the generating functions he had obtained give a finite set of 
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exact solutions of Heun's differential equation. The complete list of Carlitz results were 



gathered in JL4J], and an integral transformation connecting Heun functions with different 
parameters was derived. This last tool, combined with Carlitz results, led to an enlarged 
set of Heun functions. 

More recently we have obtained a generating function for the associated Stieltjes-Carlitz 
polynomials with rates 

A n = k 2 (2n + 2c + l) 2 , /i n = 4(n + c) 2 + fx5 n0 (3) 

A n = (2n + 2c + l) 2 , fin = 4k 2 (n + c) 2 + k 2 /j5 n0 (4) 

This result, combined with asymptotic analysis, has led to the Stieltjes transform of the 
associated Stieltjes-Carlitz polynomials |l7j 

Our aim is to show that these results lead to a finite set of exact solutions of what 
could be called the "associated Heun" differential equation which emerges as an equation 
satisfied by the generating functions of the polynomials F n (x) with the rates (|3|), Q). 

This link with orthogonal polynomials is even more fruitful since it gives a convenient 
tool to derive two new integral transformations relating associated Heun functions with 
different parameters. 

The plan of this communication is the following. 

Section 2 is devoted to a short summary on Heun's differential equation and its relation 
with orthogonal polynomials. 

In section 3 we present the associated Heun differential equation and derive two integral 
connection relations for its solutions. 

In section 4 we give a finite set of exact solutions of the associated Heun differential 
equation which generalize Carlitz results to non- vanishing (c, y). 



In section 5, in order to cross-check the Stieltjes transform obtained in JT7|], we present 
a completely different derivation which uses Karlin and MacGregor representation theorem 
and the link between orthogonal polynomials and birth and death processes. Solving the 
Kolmogorov equations leads to a Stieltjes transform which is in perfect agreement with 
asymptotic analysis. 

2 Heun differential equation 

This equation is the most general second order differential equation with four regular 
singular points 

w = 0, 1, l/k 2 , +oo < k 2 < 1 

and is described by the array of parameters 

P = {a,/?; 7, 5, e; s} a + (3 = 7 + 5 + e - 1. 

The accessory parameter s is unconstrained. We shall follow the standard notations of 



19| , p.576],||, p. 57-62] for Heun's differential equation 

w(l - w)(l - k 2 w)D 2 F + [7(1 - w)(l - k 2 w) - 5w (1 - k 2 w) - ek 2 w(l - w)]DF 

+ (a(3k 2 w + s)F = (5) 
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d 

with D = ——. 

We shall denote by Hn(P,w) the unique solution of © which is analytic for \w\ < 1 
and is normalized according to 

Hn(P,w = 0) = 1. (6) 

This equation, for arbitrary values of s, can be solved in terms of hypergeometric 
functions only for two values of the parameter k 2 : 
1) if k 2 = the solution analytic around w = is 



F 1 [ T+,T - -w) r± = a±Va 2 + s 



7 + 5-1 



7 J 2 

2) if k 2 = 1 the solution analytic around w = was shown in [F2 to be 

ft \r v ( r + a,r + (3 \ r~ 2 - 7 - a - /3 
(1 — w) 2^1 I ; it; I r = a + ycr — ap — s, a— 

For some particular values of the accessory parameter s Heun's functions degenerate 
into hypergeometric functions of the variable R(w), where R(w) is a polynomial of second 
degree in w. These values of s are listed in HlQf . 

In all what follows we shall not consider these particular cases. 

Since Hn(P, w) is analytic for \w\ < 1 we can consider it as the generating function of 
the polynomials F n (P, s), with variable s, such that 

Hn(P, w) = Yl F n( p > s V" \w\<l (7) 

n>0 

Relations (|5],|7|) imply routinely the three term recurrence relation for the F n 

(A n + p, n + 7„ - s - a(3k 2 )F n = fi n+1 F n+1 + A„_iF n _i, n > 

F_! =0, F = 1 

with 

A n = k 2 (n + a)(n + @), fi n = n{n + 7 — 1), j n — (1 — k 2 )5n. 

This recurrence exhibits the polynomial character of F n with respect to either the variable 
s or the more familiar x = s + a(3k 2 . 

One should observe on (Bp that only for 5 = do we have a true birth and death 
process with birth rate A n and death rate fi n (see || for an introduction). For 5 7^ we 
have killing in the sense of Karlin and Tavare f|,|| with rate j n . 

The finite set of exact solutions of @ obtained by Carlitz can be found in [14, p. 692] 
as well as as the following integral transform, quoted here for convenience. 

Let us define P' = {a' = 7, (3' = /3; 7' = a, 5' = 5 + 7 — a, e' = e + 7 — a; s}, then 
provided that Re 7 > Re a > 0, we have 

Hn(P, w) = —, — r f 1 dt f-^l - ty- a - 1 Hn(P t ; wt) (9) 



(8) 



B[a, 7 — Oij Jo 

for w in C\[l, +00 [. This relation, combined with the set of Carlitz solutions gives another 
finite set of solutions described in J14 , p. 693]. 
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3 Integral connection relations for associated Heun 
functions 



(10) 



Let us turn ourselves to the associated polynomials with the recurrence 

(A n + [L n + 7„ - s - (a + c)(P + c)k 2 + k 2 5c) F n = /j, n+1 F n+1 + A n _iF n _i, n > 

F_a = 0, F = 1 

with 

A n = A; 2 (n + c + a)(n + c + /3) 
< M n = (n + c) (n + c + 7 - 1) + fi5 n0 

^ 7n = (1 - + C) 

Two new parameters appear: c which is an association parameter, and // which is a co- 
recursivity parameter. Clearly, for vanishing (c, /1) the recurrence (0) reduces to (§). The 
constant terms added to s were chosen for notational convenience reasons. 
We define 

Hn(c,^P;w) = Y / Fn(c,^P)w n \w\ < 1 (11) 

n>0 

and from ([H]) it is easy to obtain 

w(l -w)(l- k 2 w)D 2 F + [(7 + 2c)(l - u>)(l - k 2 w) - 5w{l - A; 2 w) - ek 2 w(l - w)]DF 



+ 



(a + c)(P + c)k 2 w + 



c(c + 7 - 1) 
w 



W) 



8c 



c(c + 7 - 1) 
w 



fx (12) 



In view of its origin it is natural to call this equation the associated Heun differential 
equation. This departs from the usual terminology where associated differential equations 
refer to homogeneous extensions of a given differential equation whereas here we have an 
inhomogeneous extension. 

At any rate it reduces to Heun equation in the particular cases (c = 0, 11 = 0) and 
(c = 1 — 7, n — 0), according to the relations 



Hn(0,0,P;w) = Hn(P;w) 
Hn(l - 7, 0, P; w) = Hn(P; w) 



(13) 



with the array P = {1 — 7 + a, 1 — 7 + /3; 2 — 7, 5, e; s — (1 — 7)5}. 

Let us now derive two integral transformations relating associated Heun functions. To 
do this we switch from the F n to the G n defined by 



G = F , G n = Hx ■ ■ ■ fi n F n = (1 + c)„(7 + c) n F n n>l 

whose recurrence 

(A n + ii n + 7„, - s - (a + c)((3 + c)k 2 + k 2 5c) G n = G n+1 + X n -i^ n G n -i, 

G—\ — 0, Gq — 1 



(14) 



n > 



(15) 
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reveals that (— l) n G n is monic in the variable x = s + k 2 (a + c)(f3 + c) — k 2 8c. 

The basic technique to get an integral transform is to look for a mapping of the pa- 
rameters P which leaves invariant the recurrence (|15|). 

A first possibility is the mapping P' a 

a' = 7, /?' = P 

7' = a, 5' = 5 + 7 — a, e' = e + 7 — a 



s' = s, c' = c, y! — \l 



for which we have 
Using ([14]) gives 



G n (P' a ) = G n (P) n > 0. 



If Re 7 > Re a > —Re c we can write 



(c + 7) n 5(7 - a, a + c) Jo 

and inserting this in QT5|), multiplying each term by and summing n from zero to infinity 
gives 

Hn(c, /i, P, w) = — — f 1 dt t c+a -\l - ty- a - l Hn(c, u, P' a - wt) (17) 

5(7 — a, a + c) Jo 

valid for |iu| < 1. The term by term integration is allowed since the right hand side power 
series is absolutely and uniformly convergent for |iu| < R < 1. Analytic continuation 
extends this relation to C\[l, +00 [. Clearly for c = \i = we recover (^|) and there is 
another integral transformation P'^ obtained from P' a by the exchange of the couples (a, a') 
and 

A second possibility which leaves invariant the recurrence (|T3J) is P" with 
a" = 2 - a, p" = (3 + 1 - a 
7" = 7 + 1 - a , 5" = 5 + 1 - a, e" = e + 1 - a 
s" = s + (a — 1)(7 + 5 — a), c" — c + a — 1, //' = /i 

which leads to 

F n (P) = P-^F^. 
[C+l)n 

The corresponding integral transform follows analogously to fll7D 

Hn(c, u, P, w) = — r f 1 dt t c+a -\l - t)- a Hn(c + a - 1, u, P''- wt) 

B(l — a,c + a) Jo 

and is valid for 1 > Re a > —Re c and w G C\[l, +00 [. 

This second integral relation is a genuinely new result, since it changes the value of 
the association parameter from c to c + a — 1. For this reason it could not appear in the 
previous analyses where c = 0. Here too the interchange of the couples (a, a") and (3") 
leads to another mapping P'L 
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4 Exact solutions of the associated Heun equation 

Before giving the set of associated Heun functions which generalize Carlitz ones we shall 
explain, on the first of them, how they can be constructed. 
We first make the change of function 

G = w c F (18) 

which brings ([12]) to 

w(l -w)(l- k 2 w)D 2 G + [ 7 (1 - w){\ - k 2 w) - 8w(l - k 2 w) - ek 2 w(l - w)]DG 
+ [s + k 2 c(c + e + 7 - 1) + a[3k 2 w] G = c(c + 7 - l)w c ~ l + /jw c . 

The first exact solution will correspond to the parameters 

P = ja = 0, P=2> 7 = \ b= \ e= \i s = °- k2 c 2 ] c > 0. 
The change of variable 

Vw = sn(6; k 2 ) (19) 

(in what follows, concerning elliptic functions we stick to the notations of |[L9|| ; here sfw is 
the square root which is positive for real positive w) reduces the differential equation to 

d 2 G + AaG = 2c(2c - 1) (sn 2 6) c - 1 + 4//(sn 2 6) c = J \6) 

Its solution is 

G{6)=f 9 du Sin2 ^- U) j(u) (20) 
J 2-y/a 

and is meaningful provided that c > 1/2. In order to extend this integral representation 
to c > an integration by parts of (sn 2 w) c-1 is needed with the final result 

r0(w) 

w c Hn(c, /i, P; — k 2 c 2 ;w) = / du cos2^/a(9(w) — u)2c(sn 2 u) c ^ 1 ' 2 cnudio.u 

Jo 

+ / rfw V o V 1. ; ^ 4(c 2 + c 2 A; 2 + /i)(sn 2 u) c -2c(2c+ l)fc 2 (sn 2 n) 6+1 (21) 



The variable 6(w) is obtained through the inversion of relation ( 19D 

t'(w) 



l-t 2 )(l-A: 2 t 2 ) 



It is analytic for \w\ < 1; its analytic extension to the whole complex w plane has been 
described in |I|, p. 122] and is analytic but for the branch points w — 1, l//c 2 . 

We take for (sn 2 w) c-1 / 2 and w c one and the same principal branch in order to secure 
the analyticity of Hn(w) for |iu| < 1. 

It is lengthy, even if straightforward, to check that (^) is indeed a solution of ([12]) 
in the complex plane deprived with the points w — 1, 1/k 2 and that the normalization 
condition flBD does hold. 
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As mentioned in the previous section, this associated Heun function should reduce to 
a Heun function either if (c = /x = 0) or if (c = 1/2, /x = 0). In the first case a limiting 
procedure which makes use of 



f 

lim / du f(u) 2c(sn 2 ii 

c^O Jo 



/(o) 



gives 



lim Hn(c, /x, P; a; w) = cos 2yft9(w) + - (1 -cos2^( 



(7 



In the second limiting case, using relation (p0|) we get 



lim y/wHn c, /i, P; a 



A- 2 



in 



c->l/2 



sin 2^/a9(w) 



4/x 



, sin2 v /^(^H -u) 
du — -= SI1M 



2v^ 



For fj, = 0, using the relations (|I3|) we recover Carlitz results: 
Hn (o, -; -, -, -; s = a; ^ = cos (2 yfaQ (to)) 



A 3 1 1 1 + fc 2 N 

^ n 2' 1; 2'2'2 ;S = a -^ ;Uj 



sin2y / o : 6'(w) 



In order to get the remaining set of exact solutions one has to change ( |ISD into 



where A, /x, z/ take the values or 1/2. We get in this way seven more solutions to be listed 
below. 

.2' ' 2' 2' 2' 4 



P 



u>Vl - wHnic, a, P; w) = / cos 2 v / o : (6 l (w;) - w) 2c(sn 2 u) c ~ 1/2 dn 

JO 



+4(A;V + /x) 



sin 2^^) -u) 2 
art = sn w cnw 



(22) 



2v^ 



which reduces for (c = /x = 0) and (c = 1/2, /x = 0) to 

1 



113 1 , _ 

VI - wHn ( -, 1; -, -, -; s = a - w) = cos(2^/a9(w)) 

<c(i - ir)Hn \ s = ° ~ 1 ~\'^ w 



sin 2y/a9(w) 



12' '2'2'2' V 2 ; J 



: Vl — k 2 wHn(c, /x, P; if) = / cos2V?(#(w) — w) 2c(sn 2 -u) c cnw 

jo _ 

" e M sin2 v /o : (^(w) - u 



+4(c 2 + /x) 



(23) 



(fox 



2V^ 



sn u) c dnu 
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which reduces for (c = /i = 0) and (c = 1/2, // = 0) to 

/ X 113 A; 2 \ 

Vl - k 2 wHn ( -, 1; -, -, -; s = a - —;w\= cos(2^/a6(w)) 

33 13 1 2 \ an2y/a9(w 



w(l — k 2 w)Hn [1, s — a k 2 ;w 



\ ' 2' 2' 2' 2' 4 ' / 2^ 



3133 1 j 2 / \ 2 

1, -;-,-,—; s = cr k cH — 

'2 2 2 2 4 v 2 J 



r 6(w) 

(l-w)(l-k 2 w)Hn(c,/i,P;w) = / du cos2^(6{w) - u) 2c(snV c ~ J 2 

, /•*(««) , sin2^(^H-w) / 2 \ c 

+4/i / dtt — ' (sn u) cnudnw 

Jo 2-y/a 

which reduces for (c = /x = 0) and (c = 1/2, /i = 0) to 

/ 3 1 3 3 1 + A; 2 \ 

(1 -w)(l - k 2 w)Hn f 1, -; -, -, -;s = a ~' w ) = cos ( 2 V^0(w)) 

3 3 3 3 2 \ sin 2y/a9(w) 



w(l —w)(l — k 2 w)Hn -, 2; s = er — 1 — k ;w 



2' '2'2'2' ' J 2^ 

The eight particular cases for which either (c = ji = 0) or (c = 1/2, /z = 0) reproduce 
the results collected in the table |14], p. 692]. There are four other solutions: 



1,311 1 j2 . 1, 2 

-, 1; -,-,-; s = a k (cH — ) 

.2' '2'2'2' 4 v 2 } 

Let us define 

F(u) = 2c(2c + 1) (sn Vr 1/2 + 4/i(sn 2 w) c+1 / 2 
then we have _ 

y/wHn(c, n;w) — du —= F[u) 



n 



2v^ 



''2 2 2 2 v 2' 



w(l — w)Hn(c, fi;w) = [ du — F(u)cnu 

Jo 2y/o 

p = <J [ ^; ~ ~ -; s = er — - — k 2 (c + l) 2 
''2 2 2 2 4 v ; 



w(l — k z w)Hn{c, n;w) = / du ^(ujdn-u 

Jo 2^/a 

3 3 3 3 



^=^.2;^^|; a = a-l-fc 2 (c + l) a 



/— — 7^—; TT , , f e ^ , sm2y/a(8(w) - u) . 

'tml — — k z w)Hn(c, /z; u>) = / aw p= F(u)cnudau 

jo 2^/a 
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For these last cases the limiting cases (c = 0,/i = 0) and (c = 1/2, // = 0) do not 
give anything new. Furthermore the reader can check that all the solutions given here are 
correctly normalized at w — 0. 

The results (|22| , |2"3|) were first derived in [[L7j whilst all the other are new. Obviously 
this set of solutions can be further enlarged using the two integral transforms of section 3. 



5 Associated Stieltjes-Carlitz polynomials 

In |17j the Stieltjes transform of the orthogonality measure of the associated Stieltjes- 



Carlitz polynomials has been derived for the the first time. The main tool used in this 
work is Markov theorem and asymptotic analysis. It is of some interest to check this result 
using a completely different approach and this is the aim of this section. 

The strategy used here was already applied to the Stieltjes-Carlitz polynomials in 



TBI and provided for a new derivation of the orthogonality measures. Its generalization, 



described in [[TJJ and in led to a one parameter family of orthogonality measures for 



the indeterminate moment problem corresponding to the rates 

\ n = (An + l)(4n + 2) 2 (4n + 3), /j n = (An - l)(4n) 2 (4n + 1). 

What is basic in this approach is the connection between orthogonal polynomials and 
birth and death processes; the whole problem to get the orthogonality measure is reduced 
to the resolution of a linear partial differential equation. 

Let us first recall that the first family of associated Stieltjes-Carlitz polynomials are 
defined by the recurrence relation 

(A n + n n — x)F n — fj, n+1 F n+ i + A n _iF n _x, n > 

24 

F_! = 0, F = 1 

with the rates 

X n = k 2 (2n + 2c + l) 2 , fi n = 4(n + c) 2 + //5 n0 c>0 

In order to obtain the Stieltjes transform of their orthogonality measure we shall relate 
them to the birth and death process whose Kolmogorov equation is 

~^Pm,n(f) ^n—l'Pm,n~l(f) Hn+l'Pm,n+l(f) (^n l^n)'Pm,n(t) (25) 

'Pm,n(t) is the probability of a population n at time t provided that it was m at time t = 0. 
It is therefore positive and bounded by one. The link between (pf ) and ( f25| ) is provided 
by Karlin and McGregor representation theorem @,@ 



1 r°° 

V m ,n(t) = — / d^(x)F m (x)F n (x)e 

7T m Jo 



tx 



with 

A . . . A m _i 

7T = 1, TT m = , m=l,2, 

/^l • • • l-^m 



10 



G. Valent 



From this representation theorem it follows that a possible way to get \l/ is to compute the 
Laplace transform of 

V 00 (t) = / dt>(x) e~ tx 
Jo 

which we shall write _ 

V 00 (p) = / dt e- pt v 00 (t). 

Jo 

For Re p > this is nothing but 



Voo (p) 



d^(x) 



lo p + X 

closely related to the Stieltjes transform of the orthogonality measure since we have 

Voo(p) = S(-p). 

In this approach no recourse to asymptotic analysis is needed to get S(z): we just require 
^oo (p)- We shall describe in the following how this can be worked out just by solving 
linear partial differential equations. 

As a first step we consider the change of basis V m ,nit) —> Pmnif) such that 

p -" w = (I7^ p - (i) (26) 

Kolmogorov equation becomes 
dt 

P (C\\ - (1 + c )m r 

m ' nl j ~ (l/2 + c) m m "' 

with 

A n = A; 2 (2n + 2c + l)(2n + 2c + 2), fx n = (2n + 2c - l)(2n + 2c) 

whilst A n ,/i n are given by ((^). 

In order to solve (p7|) we introduce the generating function 



(27) 



ff m (f , = Vl - fc 2 w ]T P m „(t) «; n+c |u;| < 1. 

n>0 

The Kolmogorov equation becomes a linear partial differential equation for H m 

d t H m {t, w) = {Aw(l -w)(l- k 2 w)d 2 w + 2[(1 -w)(l- k 2 w) - w(l - k 2 w) 

-k 2 w(l - w)]d w } H m (t, w)-(fM+ 2 -^l) wWl - k 2 wP m0 (t) 

with the boundary condition 

H m {0,w) = \) + ° )r " w m+c VT^k^. 



(28) 
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It is convenient, from a notational point of view, to keep P m o(t); however this is related 
to the generating function H m (t,w) by 

Pmo{t) = limw-°H m (t,w) 

From now on we shall restrict ourselves to m = 0, and in order to simplify (]28f) we change 
the variable to w — sn 2 (9, k 2 ) which maps [0, 1] into [0, K\. Deleting the subscript m = 
in H we are led to 

d t H(t,9) = d 2 H(t,9) - f/i+ 2c(2C ~ 1) l (sn 2 9) c dn9P 00 (t) (29) 

\ sn z 9 J 

with 

H(0,6) = dn#(sn 2 #) c . 

Let us stress that since the moment problem for the associated Stieltjes-Carlitz polynomials 
is determined the measure \I/ is unique and therefore the solution of the Kolmogorov 



equations is unique 11 



In order to get it we extend H(t,9), a priori defined for 9 6 [0, to the interval 
9 G [-K, +K] by using the symmetry 9 <-> —9 of the equation and of the boundary value, 
and further to all real values of 9 . This last step is possible since (^) and the boundary 
conditions are periodic with period 2K. 

We shall use a Laplace transform in the variable t 

+0O 



r+oo 

H(t, 9) — ► H(p, 9)= dt e~ pt H(t, 

Jo 



to solve equation (p9|). We shall first examine what can be said on general grounds on 
H(p,6). 

Firstly since the V mn are probabilities we have the bounds 



which imply 



< Pmn(t) < n = 0,l--. t>0 



< H(t,9) < dn^(sn 2 9) C 2 F 1 ( + c ' Sn ^ 



for any real 9. From theorem 2.1 of |18, p. 38] it follows that H(p,9) is analytic in the 
domain Re p > uniformly for real 9. 

Secondly the small time behaviour of the transition probabilities is given by 

lim P mn (t) 3mn 

from this and theorem 1 of p. 181] we conclude to 



lim H(p, 9) = lim H(t, 9) = dn fl(sn 2 9Y 
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Thirdly H(p,0) must be periodic in 9, with period 2K, and for Poo(t) to exist it is 
necessary that 

]imH(t,6) = £>0. (30) 

Having stated the most useful properties of H (p, 9) let us now take the Laplace trans- 
form of equation (p9|) . We get 



d$H(p, 9) - P H(p, 9) = -H(0, 9) + A(9)P 00 (p) = J(p, 9) 

with 

A{9) = L + 2C(2C " l) \ dn£(sn 2 £) c . 
\ sn z 9 J 

This equation has for general solution even in 9 

H(p, 9) = / d<p e-^*J(p, <p) + (9^ -9) + C(p) cos\i(Jp9) 

2^/p Jo 

The integral over is convergent at provided that we take c > 1/2. 

The necessary condition ( |30|) implies C(p) =0 and the periodicity of H(p,9) in the 
variable 9 implies 

d<p e'^ J(p, 4>) = c> 1/2 



from which we deduce 

P 00{P) = —2K 



fo d(j> e-Vft A(<f>) 
Using the periodicity of A and H brings this ratio to 



d<j> cosh y/p(K - 0) H(0, , 



Poo(p) 



f K d<f> cosh ,/p(K-<j>) A{(j>) 
and the change of variable p = —z gives eventually the Stieltjes transform 

dm /* du cos y/z(K - u) H(0, u) 



rOO 



h z — s du cos \fz{K — u) A(u) 

so that if we define 

D(c,fx;z)= J du cos \fz{K — u) \2c(2c — 1) + fisn 2 uj dnu 



r(2c+ 1) 



we end up with 

*.,»(.) D(c+1,M 



lo z — s D(c,p,;z) 

in perfect agreement with the result derived in [|17]] . In this reference the polynomials with 

rates 

A n = (2n + 2c + l) 2 , fa = 4:k 2 (n + cf + fik 2 5 n0 < k 2 < 1 

have been shown to follow from the result (|3l|) using the transformation theory for Jacobian 
elliptic functions. 
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